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Versatile controllability of interactions in ultra¬ 
cold atomic and molecular gases has now reached 
an unprecedented era where quantum correla¬ 
tions and unconventional many-body phases can 
be studied with no corresponding analogs in solid 
state systems. Recent experiments in Rydberg 
atomic gases have achieved exquisite control over 
non-local interactions HHg> allowing novel quan¬ 
tum phases unreachable with the usual local in¬ 
teractions in atomic systems. Here, we study 
Rydberg dressed atomic fermions in a three di¬ 
mensional optical lattice predicting the existence 
of hitherto unheard-of exotic mixed topological 
density wave phases. We show that varying spa¬ 
tial range of the non-local interaction leads to a 
rich phase diagram containing various bond den¬ 
sity waves, with unexpected spontaneous time- 
reversal symmetry breaking. Quasiparticles in 
these chiral phases experience emergent gauge 
fields and form three dimensional quantum Hall 
and Weyl semimetal states. Remarkably, certain 
density waves even exhibit mixed topologies be¬ 
yond the existing topological classification. Ex¬ 
perimental signatures of density waves and their 
topological properties are predicted in time-of- 
flight measurements. 

Topological quantum states of matter have attracted 
considerable theoretical and experimental attention in 
condensed matter research in the last decade, aiming at 
robust physical properties protected (often due to some 
underlying symmetry) against variations in microscopic 
details. For non-interacting fermions, topology of gapped 
phases has been classified for arbitrary dimensions 00 - 
Examples are Z class Chern insulators in two dimensions 
manifesting quantum anomalous Hall effect and 

time-reversal invariant Z 2 topological insulators in three 
dimensions with robust magnetoelectric response HBHS]. 
For classification of gapless systems, no unified theoreti¬ 
cal framework has emerged yet, nonetheless recent stud¬ 
ies on semimetals suggest one theoretical route by char¬ 
acterizing topological properties of band touching points, 
such as flux monopoles nmg. In this paper, we provide 
one fascinating possibility beyond the present theoretical 
paradigm in an interacting system, with mixed gapped 
and gapless topologies emergent in a three dimensional 
bond density wave state. We demonstrate the existence 
of such an exotic topological density wave in a Rydberg 
dressed atomic Fermi gas via microscopic calculations. 
We also propose an experimental scheme to extract the 
topological properties based upon time-of-flight signals. 

We consider Rydberg dressed atomic Fermi gas in 


a three dimensional (3d) cubic optical lattice. This 
Rydberg atomic system has density-density interac¬ 
tions described by a non-local potential, V (r) = 

V 6 /[l + (|r|/r c ) 6 ], where V§ describes the interaction 
strength, and r c the interaction range determined by 
the Condon radius in Rydberg dressing |Ti7112U] . Con¬ 
sidering a deep lattice, the kinetic motion of atoms arises 
mainly from quantum tunneling between nearest neigh¬ 
boring sites providing a unique kinetic energy scale, t. 
We will focus on a lattice with octahedral Oh symmetry. 
The intrinsic properties of this system will depend on 
only two dimensionless parameters— V^/t and r c /a (with 
a the lattice spacing), representing the strength of inter¬ 
action and its non-locality, respectively. 

At half filling, namely one atom per two lattice sites, 
the Fermi surface is perfectly nested in the deep lattice 
limit, because the single particle energy dispersion sat¬ 
isfies £k+Q = — £k, with Q = (7r,7r,7r) (see Methods). 
The susceptibility towards forming particle-hole pairing 
or density waves 

p k = (^(k + QMk)), 

(with 0(k) a Fourier transformed annihilation operator) 
is logarithmically divergent due to nesting and the sys¬ 
tem has generic instabilities with repulsive interactions. 
With nearest neighbor interactions for spinless fermions, 
the ground state density wave order pk is momentum 
independent, implying a trivial 3d checkerboard pattern 
in real space. However, for Rydberg dressed atoms, the 
interaction range, r c can be comparable to several lat¬ 
tice spacings, which then frustrates the checkerboard pat¬ 
tern, giving rise to possibilities of unconventional density 
waves [2T1122] . 

Within one-loop renormalization group (RG) analysis, 
we find that the strength of density wave instability is 
determined by the eigenvalue problem of a rescaled in¬ 
teraction matrix 7 whose symmetry is Oh xT (see Meth¬ 
ods). Different density waves correspond to eigenvectors 
of 7, which are classified according to irreducible repre¬ 
sentations of the symmetry group [23] . The symmetry 
classification of density waves and their representative 
irreducible basis functions are shown in Supplementary 
Table 1. The trivial density wave is A^~ g and all other 
density waves necessarily have more complicated momen¬ 
tum dependence for symmetry reasons. The latter leads 
to particle-hole pairing on the bonds in real space Ell- 

Before presenting our microscopic results, let us first 
look at a nontrivial density wave, Tf~ u +T^, which has a 
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FIG. 1. Density wave topology. The color index the phase angle of the momentum dependent density wave order Ak projected 
on the Fermi surface, b and c use the the same color scheme as shown in a. a, b and c correspond to topologically distinct 
density wave states. The red dotted lines are vortex lines around which the phase of Ak changes by 2mn (m 0). Right 
on vortex lines, Ak vanishes. In a, one vortex line crosses the Fermi surface. In both a and b, all vortex lines cannot be 
adiabatically removed without touching the Fermi surface, and thus are topologically robust, whereas in c, the vortex line is 
adiabatically removable and is thus topologically trivial. 


complex order Ak (see Methods) of the form, 

Ak = A t + [sin k x (cos k y — cos k z ) + sin k y ( cos k z — cos k x )} 

+iV 2A t - sin k z . 

1 1« 

Here A r - and A T + take real values. Such a density 
wave order has a rich topological structure in 3d momen¬ 
tum space (Fig. [TJi), with three vortex lines located at 
(0, Z, 0), (Z, 0, 0) and (1,1, 0), with l £ [— n/a, 7r/a). These 
vortex lines are topologically robust in the sense that 
they cannot be smoothly removed without touching the 
Fermi surface. The Bogoliubov-de Gennes Hamiltonian 
to describe quasiparticles in the density wave background 
reads, 


-ZZBdG (k) 


£k A k 
K -Ck 


( 1 ) 


which can be rewritten in terms of Pauli matri¬ 
ces as i?BdG(k) = h(k) • a. At the points kj = 
(±27r/3, ±27 t/ 3, 0) where vortex lines cross the Fermi 
surface, the magnitude of h vanishes, leading to zero 
energy quasiparticles. Around each nodal point, say 
k = k,| + (6k x , Sky, 8k z ), we have h x ss |A T + ( Sk x — 5k y ), 

h y ps — \/2A T - <5fcs and h z ps y /3 t(Sk x + 5k y ). The 
quasiparticles near the gapless points are thus Weyl 
fermions HI US] with highly anisotropic velocities. Be¬ 
sides this gapless topological property, it is worth noting 
that the other two vortex lines that do not cross the Fermi 
surface are also topologically robust. 

In general, we argue that the topology of 3d density 
waves can be characterized by two integers ( 711 , 712 ), 711 
being the number of vortex lines that cross the Fermi sur¬ 
face, and 77.2 the number of other vortex lines that cannot 
be contracted to one point without touching the Fermi 
surface. The topological numbers ni and 77.2 thus encode 
gapless and gapped topologies, respectively. We empha¬ 
size that any two vortex lines connectible by the Q vector 


are equivalent. For gapped density waves, n\ is necessar¬ 
ily zero, and we have a Z classification, whereas for gap¬ 
less states, both 77i and 712 can be finite and the classifica¬ 
tion is Z x Z. Using this scheme, the T^ u + T£ u state has 
topological numbers (1,2), exhibiting mixed topologies. 
The density wave topology as constructed is measurable 
by time-of-flight techniques in atomic experiments (see 
Methods). 




FIG. 2. Leading density wave instabilities for Rydberg 
dressed atomic Fermi gas. a, most dominant eigenvalues of 
the 7 matrix (see Methods), representing the strengths of 
particle-hole paring instability of corresponding channels. In¬ 
creasing 77, the A+ g channel gets suppressed, and other chan¬ 
nels, Tj“, E~, T^ g , and T% u , with nontrivial momentum de¬ 
pendence are enhanced, b, step-like interaction form for Ry¬ 
dberg dressed atoms. Two atoms far apart behave as in Ryd¬ 
berg excited states and interact with van-der waals potential, 
whereas they become more ground state like at short distance, 
c, the Fourier transformed interaction V'(Q) [Q = (77,77,77)]. 
As we increase the cut-off r max (see Methods), V’(Q) con¬ 
verges and the truncation error is already negligible when 
7*max/7 'c ~ 3. 
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We now discuss the actual density waves supported 
by Rydberg dressed atoms. With infinitesimal repulsive 
interaction, the relative strengths of instabilities in dif¬ 
ferent channels are fully determined by the non-locality 
strength r c /a. The leading instabilities, as quantified by 
the corresponding eigenvalues of the 7 matrix, are shown 
in Fig. [2] When r c is small, the interaction is essen¬ 
tially nearest neighbor, for which the 7 matrix is mainly 
determined by the momentum independent part of inter¬ 
action V"(Q) (the Fourier transform of the interaction at 
Q). This makes the trivial Af g density wave completely 
dominant over all other density wave states. Increas¬ 
ing r c , the longer-ranged part of the interaction becomes 
more important, which suppresses V'(Q) (see Fig. [ 2 J 3 ) 
and consequently the Af g density wave. Simultaneously, 
other density wave channels with non-trivial momentum 
dependence get enhanced. The dominant instabilities are 
Ti“ , E~ , T^ g and T ^ u , and we have a transition from Af g 
to Ti u around r c /a = 1.46. The phase diagram based 
on the instability analysis from RG is suggestive, and is 
rigorous in leading order. However, considering finite in¬ 
teractions, non-linearity may lead to significant physical 
effects. In particular, when r c reaches the scale of lattice 
spacing, the instability strengths in different channels are 
quasi-degenerate, which we attribute to the “step-like” 
feature of the Rydberg dressed interaction (see Supple¬ 
mentary Information). Non-linear effects must then be 
taken into account to determine the actual phase dia¬ 
gram. 

To incorporate nonlinearity, we numerically simulate 
a system of size 64 x 64 x 64 with periodic boundary 
condition using self-consistent methods (see Methods). 
The solution for the ground state density wave order Ak 
is expanded in terms of basis functions 


In addition, we find three gapped phases, A+ g , A+ g + T lu 
and E~ +T^ g , whose quasiparticle density of states has a 
“hard-gap” (Fig. [3]b). In atomic experiments, the quasi¬ 
particle spectrum and density of states can be directly 
measured by radio-frequency spectroscopy [241 [25] . 



FIG. 3. Phase diagram for Rydberg dressed atomic Fermi gas. 
In this plot, the interaction strength is fixed to be V%jt = 4. 
To quantify contributions from different classes, we introduce 
A q = IAo:,/?! 2 , (see equation §). a shows the phase 

diagram with varying Condon radius r c . As we increase r c , 
we get a sequence of density wave phases, A+ g , A+ g + 7j“, 

T 2g + E g ! E g + T lu + T 2g + T 2u and T lu + T 2W T l le firSt 

three are fully gapped and the last two are gapless. The 
phase transition between the two gapless phases is found to be 
second order. The corresponding broken symmetry is space 
inversion. The phase transitions among the gapped phases 
are first order, b and c show the density of states for gapped 
and gapless phases, respectively. 


Ak = ^ A Qi/ 3 ^ a ^(k), ( 2 ) 

o:,/3 

where (t> a ,p( k) are the basis functions from symmetry 
classification (see Supplementary Table 1), with a label¬ 
ing different classes and /? different functions within one 
class. In the self-consistent theory including non-linear 
effects, spontaneous symmetry breaking could occur and 
superpositions of density waves from different classes are 
allowed, subject to the constraint Ak = A£. + q. The co¬ 
efficients A a> 0 are required to be real. 

In our self-consistent calculations, spontaneous sym¬ 
metry breaking of Oh x T is found for a large range of r c 
and the symmetry-broken ground states support various 
superposed density waves, yielding a rich phase diagram 
(Fig# Surprisingly all the superposed states are found 
to break T symmetry and consequently the time-reversal 
symmetry due to the resultant complex structure in Ak- 
The topological gapless density wave T]” + T^ u is indeed 
an allowed solution as argued earlier. The other gapless 
phase that occupies a large region in the phase diagram is 
E~ + X)” + T,f g + T 2 + . In both gapless phases the quasipar¬ 
ticle density of states has a “soft-gap” feature (Fig. § 0 . 
a signature of the emergent topological Weyl fermions. 


We now discuss topological properties of the density 
waves in the phase diagram. In E~ + T^ u + T^g + T£ u 
state, as compared to T± u + T^~ u , the density wave 
gets additional contributions. The vortex lines are no 
longer straight, but otherwise have the same topology as 
shown in Fig. These two gapless density waves are 
then topologically equivalent, having the same topolog¬ 
ical numbers (1,2) in our classification scheme. In the 
E~ + T^ g state which is gapped, the vortex line (with 
vorticity 2) is located along one axis (Fig [l}o) , assumed 
to be the k z axis in spontaneous symmetry breaking. 
This state has topological numbers (0,2). Actually in 
this state, the quasiparticle Hamiltonian Hq dG(k), with 
arbitrary fixed k z , has finite Chern number 2. Such 
a state thus exhibits 3d quantum anomalous Hall ef¬ 
fect, featuring chiral surface modes, which mediate highly 
anisotropic transport properties—dissipationless in x or 
y direction but diffusive in the z direction mum,. In the 
other gapped state, A+ g +T-^ U . the vortex line is shown in 
Fig.jlJ.’, but such a vortex line is unstable. The A\ g +T-^ U 
state thus has trivial topological numbers ( 0 , 0 ), topolog¬ 
ically equivalent to Af g . 

We expect that our proposed topological density waves 
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to be generic for Rydberg dressed atomic Fermi gases 
even if the interaction potential is not precisely of the 
specific r 6 form. The key qualitative ingredient is the 
step-like feature in the interaction. To confirm this, we 
also have carried out calculations for Rydberg p -wave 
dressed atoms [28J, and similar topological states are in¬ 
deed found (see Supplementary Information). We there¬ 
fore believe that our prediction for novel collective topo¬ 
logical density waves should apply generically to dressed 
Rydberg optical lattices. 

Methods 

Density wave instability under infinitesimal re¬ 
pulsion. The single-particle energy dispersion for 
the 3d cubic lattice with Oh symmetry is e k = 
—2 1 (cos k x + cos k y + cos k z ). Having non-local density- 
density interaction with potential V (r—r'), the scattering 
vertex among single-particle modes ki, k 2 , k 3 and k 4 is 

T(k 1; k 2 , k 3 , k 4 ) = \ {W(k 2 - k 3 ) + R(ki - k 4 ) 

-H(k 2 -k 4 )-H(k 1 -k 3 )}, 

with U(q) the Fourier transform of U(r). At half fill¬ 
ing, from the Fermi surface nesting, repulsive interactions 
cause instabilities in density wave channels. Such insta¬ 
bilities can be tracked by renormalization group (RG) 
flow of the effective couplings near the Fermi surface, 


wave order are given by 

1 


£ k = 


Ak = 


U(k-q+Q)-U(k-q) 
Pq. 


V(Q)-R(k-q) 

q 3 J 

with Jq = f 4$’ n q = (' i /' t (q)V’(q)) i and p q = 
(ipUq + Q)ift(q)) . The self-consistent equations are ob¬ 
tained to be 


Ak = 


1 


V(k~q)~V(Q) A q /|e q |, 


Sk — 2 / ^( k q) [e q + £q]/|e q |, 


with quasiparticle spectra e q = ± y/(e q + £ q ) 2 + |A q | 2 . 
With moderate interaction strength, self-energy correc¬ 
tions could be significant, but the modification to the 
Fermi surface is reasonably weak (see Supplementary Fig¬ 
ure 1). In numerics, we introduce a long-range cutoff 
r ma j, so that U(|r| > r max ) is set to be zero. For all nu¬ 
merical results presented in the main text, r max is set to 
be 6 x a and varied to ensure that results do not change 
qualitatively. 


Experimental signature of density wave topology. 

In experiments, the momentum dependence of the den¬ 
sity wave order can be extracted from time-of-flight mea¬ 
surements. Using mean field theory, it is straightforward 
to show that 


9kf,k' f = r(k'. + Q, k f, k' f , kf + Q), 

The RG equation at one-loop level [29, ,30] reads 

- AdA7 k/ .k^. = j ( 27 T ) 2 ^kf .q/ 7q/ ,k' f ’ 

with 7 a rescaled interaction matrix, 

lk f y f = [^Vf(k f )]- 1/2 gk f y f [2nv f (k' f )\ 1/2 , 

where Vf(kf) is the Fermi velocity. In the RG flow, the 
largest eigenvalue of the 7 matrix diverges most quickly 
as approaching the low energy limit. The magnitudes of 
eigenvalues quantify the relative strength of instabilities 
in different channels. The symmetry group of 7 is Oh x 7~, 
with T a two-element group involving the transformation 
T(k)=k + Q. 

Non-linear effects and self-consistent theory. The 

variational state we choose to describe density wave order 
is [21 

W = Ilk [wkV ,t ( k ) +v k ipi{k+ Q)] |fi), 

with |fl) the vacuum, and k running over one half of the 
Brillouin zone, e.g., k z > 0. The self-energy and density 


|pk | 2 = (n k )(nk+Q) - (n k n k +Q). (3) 

At finite temperature, to extract |pk| 2 j it is necessary to 
measure both the product (n k )(n k +Q) and the correla¬ 
tion (rikUk+Q). At half filling, the latter originates purely 
from thermal fluctuations. The former is dominant when 
the temperature is not too high as compared to the quasi¬ 
particle energy gap (or the “soft-gap” for gapless states). 
This implies that it should be experimentally straight¬ 
forward to get a precise measurement of |pk| 2 - I n the 
zero temperature limit, the density wave order is com¬ 
pletely determined by the product |pk | 2 = (n k )(n k +Q). 
Moreover from the relation, 

Pk = bkF(£k) - n F (£k+Q)] (4) 

l e k 

the vortex lines of Ak immediately follow by tracking the 
vanishing points of pk | • With vortex lines located in the 
3d momentum space, the topological numbers of density 
waves can be easily extracted. 
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SUPPLEMENTARY INFORMATION 

Exotic topological density waves in cold atomic Rydberg fermions 

S-l. A CONSISTENT THEORY FOR DENSITY WAVE ORDER 

The model Hamiltonian to describe a Rydberg dressed atomic Fermi gas in a three dimensional (3d) cubic optical 
lattice is H = H 0 + H lnt with 


Hq ^ ' [ la^rb+ae Q T h.C. 




f f 

c'c , c v >c r . 


(Si) 


r,r 7 


Here a is the lattice constant, a index the three directions x, y and z, and c r the fermionic annihilation operator 
positioned at r. The atomic interaction potential, assuming the Rydberg state is s-wave, takes an isotropic form, 
V{v) = V§/[1 + (|r|/r c ) 6 ], for r^O. For fermionic atoms considered here, F(0) is set to be zero. 

To describe a density wave ground state, we take a variational state 


l*>=IIk 


Wk4 + fk4 +Q 


|H). 


(S2) 


Here |fi) is the vacuum, and we have a normalization condition |w k | 2 + |fk | 2 = 1. In fL, k runs over one half of the 
Brillouin zone, say k z > 0. We can define u k and u k in the whole Brillouin zone by the relation iik+Q = Vk, t’k+Q = u k. 
Half filling is already enforced within this variational ansatz. The variational energy cost E = ('k|i?o + -ffintl’I') is 
obtained to be 

1 r ,_ 

r\/(ck + S k ) 2 + |A k |2 


E = — 


(2t r) 3 
d 3 k d 3 k' 


1 

+ 2 J (2t r) 3 (2 t r) 3 


A k EQ 1 (k — k')A k / 


(S3) 


with Vcj(k — k') = F(k — k') — F(Q). Minimizing the energy, we get the Bogoliubov-de Gennes (BdG) equation 


with 


where the self-energy is 


#BdG(k) 

I^BdG(k) = 

1 f d 3 q r 


«k 

t’k 


= £k 


Uk 

Vk 


e k + S k 

AT 


Ck — Si 


2 J (2ir) 3 

and the off-diagonal term is the density wave order, 

d 3 q 


y(k-q+Q)-V(k-q) 


A k = 


(27t) 3 


F(Q)-Uk-q) 


(S4) 


(S5) 


(S6) 


(S7) 


The chemical potential has been shifted by |F(0) to compensate Hartree corrections to the self-energy. By definition, 
the density wave order satisfies A k = A£ + q, which is an important difference of particle-hole pairing from particle- 
particle paring in superconducting states. Diagonalizing the BdG Hamiltonian, quasiparticle energies are given as 
£ k = ±y/(e k + S k ) 2 + |A k | 2 . In variational calculations, we choose the negative branch, and the self-consistent 
equations are then obtained to be 

A k = \ ( [y(k-q)-E(Q)l 


E k = - 
2 


(2tt) 3 
d 3 q 
(2tt) : 


Aq/|£q|, 


F(k-q) [e q + E q ]/|e q |. 


(S8) 
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The effective potential for A k , defined by C/[A k ] = (2£[A k ] — E[0]) /N a , reads 


U[ A k ] = -i 




d 3 k 
(2t r) 3 
d 3 k 


Ck 


d 3 k' 


(27t) 3 (27t) 


a tvr 1 


H k ) 2 + |A k | 2 — 
q (k — k')A k /. 


e k 


S k |} 


(S9) 


With weak interaction, the density wave order A k is restricted to the region near the Fermi surface. We thus rewrite 
the momentum in terms of Fermi momentum and the distance from the Fermi surface as 


k = k f + le±, 


(S10) 


where e± is a unit vector perpendicular to the Fermi surface. The effective potential now reads 

721 


U[ A k ] = -\ 


f d 2 kf 

[ A dl \ 

j (27t) 2 

Jo 2tt i 


■v 2 f (k f )l 2 + \A kf \ 2 - v f {k f )l 


d 2 k f d 2 k! f ^ 

(2tt) 2 (2tt) 2 k /' 9k / k '/ k /’ 


(Sll) 


where the coupling matrix [ g\ is the matrix [Vq] projected onto the Fermi surface. To see the momentum dependence 
of the A k , we expand t/[A k ] to the second order as, 




d% |A k ,| 2 
(2tt) 2 v f (kf) 


f d 2 kf d 2 k' f 

l- 

•+-, 

<1 

1_ 

J (2t r) 2 (2?r) 2 

.vVf( k /)_ 


^/u/( k /)5 k/ 1 k^ t '/( k /) 


Aw/ 


v f( k f 


,(S12) 


where an infrared cutoff r is introduced to regularize the perturbative expansion. This scale r can be physically 

identified to be temperature. Minimizing the effective potential, , k/ is an eigenvector of the matrix 

V v f( k f) 


\/ r fikf)tj k : k ' f y/vf( k' f ) = 7 kfk - (S13) 

with the minimal eigenvalue, or equivalently the eigenvector of 7 kfkf ' with the maximal eigenvalue A ma x- The corre¬ 
sponding transition temperature is estimated to be 

t c = Aexp(—27rA“ ax ), (S14) 


which is valid for infinitesimal interactions. 

It is worth noting that the couplings . 9 kjk ' actually depend on the cutoff A. The effective potential U[A] is, on the 
other hand, independent of A [30], i.e., 


Ad A C/[A k ] = 0, 

from which the renormalization group /3-function follows 

/(_j2 

^^5k /q/ [27n; / (q / )]- 1 5 q/k ^. (S15) 

This renormalization equation can also be derived by momentum shell scheme [29] . Analyzing the renormalization 
group flow, the momentum dependence of the density wave order can be determined from leading divergent channels 
in the renormalization group flow, analogous to the procedure of extracting dominant channels in unconventional 
superconductivity m- This approach is equivalent to minimizing the effective potential at quadratic order. 


S-2. AN APPROXIMATE CRITERION FOR UNCONVENTIONAL DENSITY WAVES TO DOMINATE 


As discussed in Section |S-1| the leading density wave instability is determined by the eigenvalue problem of the 7 


matrix (equation (S13)), which in general requires numerics. Here we would like to derive an approximate criterion 
for unconventional density waves to occur, and give some intuition why such orders are supported by Rydberg dressed 
interactions, assuming that the momentum dependence of fermi velocity is negligible. Under this assumption, the 
eigenvalue problem of 7 becomes equivalent to that of 


V(k —k')-U(Q). 
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Supplementary Figure 1. Illustration of self-energy corrections. In this plot we choose r c /a = 1.86 and Ve/t — 4. a, b 
and c show contour plots of bare dispersion ek with k z a = 0, \ and , respectively, d, e and f show tk + £k f° r the same set 
of k z . Although self-energy corrections are significant, its modification to the Fermi surface (defined by £k + £k = 0) i s weak. 
The property holds for other choice of parameters as well. 


From Fourier transformation, we know 

r d 3 k' 

J j2ny> [V{k k ' } " ' r = - ^(Q)^,o]e- ik ' r . 


(S16) 


The density waves, Ak, are then given by the plane waves e _,kr , r labeling different solutions. The eigenvalue 
associated with the trivial density wave with r = 0 is then — V(Q). For other density waves having momentum 
dependence, with r ^ 0, the eigenvalues are V(r). With isotropic interactions, the density wave solutions with the 
same |r| are degenerate. In particular with step-like interactions as for Rydberg dressed atoms, the solutions with 
|r| < r c are actually all quasi-degenerate. This would strongly amplify non-linear effects, giving rise to possibilities 
for superposed density waves. 

The criterion for unconventional density wave to dominate over the trivial one is, 

max{R(r)} > -R(Q), (S17) 


which implies for the Rydberg interaction that 


1 

1 + (a/r c ) 6 


>-£ 

r/0 


(_1 Y^/ a + r y/ a +r z /a 

i + (MAc) 6 ' 


(S18) 


From Supplementary Figure [ 2 } the parameter range r c /a £ (1.45,2.0) should be ideal to look for unconventional 
density waves, which is precisely the region we focus on in the main text. For general interactions, the most likely 
region to look for unconventional density waves can be easily identified by the derived criterion in equation (S17). 
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r c /a 

Supplementary Figure 2. Estimation of relative instability strength towards trivial and non-trivial density waves. The 

solid and dashed lines correspond to 1+ ^/ r )« an< i — ; i+( | r |/ r ) 6 -> respectively, representing instability strengths 

of momentum independent and dependent density waves. The latter is significantly dominant in the region, r c /a € (1.45,2.0). 


S-3. NUMERICAL TRICKS IN SOLVING SELF-CONSISTENT EQUATIONS 


In self-consistent iteration (equation (S8|), a brute force treatment would lead to numerical cost of 
the number of lattice sites, which makes the numerics too heavy for large systems, say with N s = 64 3 . We 
the problem by rewriting 


with N s 
simplify 


A k =i^e- lkr U(r) 

r 

£k = ^e-‘ kr V(r) 


d 3 q 
(2t r) 3 

d 3 q 

( 2^)3 


r A q /|e q | 


~nV(Q) 


e iq r (e q + S q ) /|e q | 


d 3 q 

(2tt) 3 


A q /|e q |, 


(S19) 


Following the bracketing order above in numerical iteration, the cost is greatly reduced to O ((T’ max /a) 3 IV s ), (?' max is 
the long-range cutoff, fixed to be 6 x a in calculations), which makes it feasible to simulate large systems. 


S-4. SYMMETRY CLASSIFICATION OF DENSITY WAVES 

In this section we give the classification of density waves according to Oh * T (see Methods in the main text). 
The basis functions for density waves with different symmetries can be constructed by a standard projection method 
in group theory j23j . Classification of density waves and representative irreducible basis functions are shown in 
Supplementary Table □ The basis functions are orthonormal. Respecting all symmetries, superpositions of density 
waves from different symmetry classes are not allowed. 


S-5. TOPOLOGICAL PROPERTIES OF DENSITY WAVES 

Details of topological properties of density waves in the self-consistent phase diagram (Fig. 3 in the main text) are 
provided in this section. In E~ + T^ a + r T,^ g + T^ u state, as compared to T± u + T^, Ak gets an additional contribution— 
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4 + 

9 

A ia 

A tu 

A Tu 

1 , cos k x cos k y + cos k y cos k z + cos k z cos k x , 

i (cos k x + cos k y + cos k z ) 




A 2q 

Atu 

A 2 u 





K 

K 

E+ 

Eu 

J cos k z (cos k x — cos k v ) 

1 2 cos k x cos k y — cos k z (cos k x + cos k y ) 

J i (cos — cosfcj,) 

1 i (2 cos k z — cos k x — cos k y ) 



T+ 

1 ig 

T u 

T+ 

1 1 u 

Tr u 



\ 

sin k x (cos k y + cos k z ) 
sin k y (cos k z + cos k x ) 
sin k z (cos k x + cos k y ) 

I i sin k x 
\ i sin k y 
[ isinfc. 

T+ 

1 2 g 

T 2q 

T + 

± 2u 

r ^2u 

[ sin k x sin k y 
s sin k y sin k z 
( sin k z sin k x 


1 

sin k x (cos k y — cos k z ) 
sin k y (cos k z — cos k x ) 
sin k z (cos k x — cos k y ) 



Supplementary Table 1 . Classification of three dimensional density wave orders according to irreducible representation of 
the symmetry group Oh x T. In the labeling of different classes, the superscript ± tells how the density wave transforms under 
T symmetry. The basis functions representing particle-hole paring of upto next nearest neighboring sites are given. Due to the 
constraint Ak = Ak + q, the basis functions in T even (odd) classes are real (purely imaginary). The normalization constants 
for the basis functions are neglected in this table to save writing. 


y/2A T + sin fc z (sin k x — sinfcj,) + iA„- [2 cos k z — cos k x — cos k y \ /-\/3. Treating these additional components, 
A e -/A t - and A T + /A T + as perturbations, we get three vortex lines given by ^0, l, A E - / (y/6A T - ) (cos l — 1)^ , 
(/,0,A B -/(v^A r - ) (cos l — 1)^ , and {l,l,k z = 2A E -/(v / 6A T - ) (cos? — 1)^ , the third of which crosses Fermi sur¬ 
face. The E~ + 7)“ + T^g + T 2u state is topologically equivalent to Tu, + Tl^, having the same topological numbers 
(1,2). 

In the E g + T 2g state which is gapped, Ak takes a form 


Ak « 2A t + sin k x sin k y + iA E - (cos k x — cos k y ). 

The resulting vortex line (with vorticity 2) is located along the k z axis (note that the Oh symmetry has been 
spontaneously broken), thus not crossing the Fermi surface (see Fig. [Tju in the main text). This state has topological 
numbers (0, 2) in our classification scheme. In the other gapped state, Af g + T± u , we have 

A k « A a + i - 2/v / 3A a + 2 ^ cos k a cos k a ' 

+ iy/ 2 A t - sinfc z . 

The vortex line of this state is shown in Fig. [TJ: in the main text, but it is unstable because it can be adiabatically 
contracted to one point without touching the Fermi surface. The A+ + T)” state has trivial topological numbers 
( 0 , 0 ), topologically equivalent to A^ g . 


S-6. DENSITY WAVES WITH RYDBERG p-WAVE DRESSED ATOMS 

In this section we present results for density waves in p -wave dressed Rydberg atoms |28| , whose interaction takes 
an anisotropic form 

V(r ± 0) = ^ sin 4 (0), (S20) 

1 + (|r|/r c ) b 

with 6 the polar angle. The symmetry group to classify density waves in this system is then D 4^ x T. The symmetry 
classification is shown in Supplementary Table 0 
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A + 
g 

A ia 

A tu 


1 , cos k x cos k y , (cos k x + cos k y ) cos k z 

i( cos k x + cos k v ), i cos k z 

(cos k x + cos k y ) sin k z 

i sin k z 

Et 

Ey 

E+ 

E u 

I sin k x sin k z 


I sin k x cos k y , sinfc^cosfcj 

I i sin k x 

1 sin k y sin k z 


1 cos k x sin k y , sinfcycosfc- 

1 i sin k v 

K 

B u 

EL 

B lu 

(cos k x — cos ky) cos k z 

f (cos k x — cos k y ) 

(cos k x — cos k y ) sin k z 


K 


to 

B 2u 

sin k x sin k y 





Supplementary Table 2 . Classification of density wave orders for non-local interactions up to next nearest neighbor according 
to irreducible representation of the symmetry group x T. 


For simplicity, we truncate the long-range part of the interaction by taking r max = \/2 x a. The interaction is then 
modeled to be 


V(±x) = V(±y) = V n „ 

V(±x ± y) = V^ nn 

V(±x ±z) = V(±y ±z) = V: m . (S21) 


All other terms are taken to vanish. With tunability of the lattice g eome tries, the ratio F r f nll /K m and V^ m JV rm is 
largely tunable between 0 and 1. The approximate criterion (equation ( S17|) ) for unconventional density waves in this 
model becomes 


3F„„ < 4V^ nn + 8F n z nn . (S22) 

We then numerically solve self-consistent equations for p-wave Rydberg dressed atoms. Both symmetric (t z = t x = t y ) 
and asymmetric ( t z = 4 t x = 4 t y ) tunneling cases are studied. 

a b 




Supplementary Figure 3. Phase diagram from self-consistent calculations with for p -wave Rydberg dressed atoms. In this 
plot we choose V nn /t z = 2, Kfnn/Fnn = 0.2. a and b show the symmetric ( t z = t x = t y ) and asymmetric ( t z = 4t x = 4t z ) 
tunneling cases, respectively. The A]~ g component A a + is rescaled by a factor 4 for illustration purposes. 


For the symmetric case (see Supplementary Figure |3jr) , we find Af g , Af g + B lg , B lg and B lg + B^g density waves. 
The Fermi surface is fully gapped out for both Af g and Af + B^ g states. In the BA state the Fermi surface is only 
partially gapped out, and we have gapless lines (protected by T symmetry), from which the density of states at low 
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energy has linear behavior, i.e., D(e) oc e. In the B± + B^g state, the Fermi surface for this state is fully gapped out, 
and the complex order A k has a vortex line (with vorticity 2) elongated along the k z axis. This B~[ g + B^g state is 
topologically equivalent to the discussed E~ + T^ g state in the main text. 

For the asymmetric case (see Supplementary Figure [3j:>), we find A+ g , A^ g + B^ g , A^ g + B^ g + E+ density waves. 
In the A^ + B^ g state, T is not broken and we have gapless lines. In the A^[ + B^ g + E+ state, the density wave 
order takes a form, 


A k = iA 4 - (cos k x + cos k v ) + iA B - (cos k x — cos k y ) + A E + sin k x cos k. 


(S23) 


which has two straight vortex lines along the k z direction, given by (k x = 7r/2, k y = ±7r/2). These vortex lines cross 
the Fermi surface, and we have four Weyl nodes. Topological numbers for this state are (2,0). 

For both cases, the transition from trivial state to other density waves roughly occurs at V^ nn /V nn = 0.35, with 
Vnnn/Vnn fixed at 0.2, which agrees with the estimate from equation (S22). 



